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Abstract 

We discuss the results obtained by fitting the lattice data of the gauge-invariant field 
strength correlators in QCD with some particular functions which are commonly used 
in the literature in some phenomenological approaches to high-energy hadron-hadron 
scattering. A comparison is done with the results obtained in the original fits to the 
lattice data. 
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1. Introduction 



The gauge-invariant two-point correlators of the field strengths in the QCD vacuum are 
defined as 

V mva (x) = g 2 (0\Tr{G w (0)S(0,x)G U(T (x)S^O,x)}\0) , (1.1) 

where G w = T a G^ p and T a are the matrices of the algebra of the colour group SU(iV c ) 
in the fundamental representation (For N c = 3, T a = A a /2, where A a are the Gell-Mann 
matrices). The trace in (1.1) is taken with respect of the colous indices. Moreover, in Eq. 
(1-1), 

S{0,x) = Pe^p(ig£dz' x A fl {z)j , (1.2) 

with Ap = T a A^, is the Schwinger phase operator needed to parallel-transport the tensor 
G ucr (x) to the point 0. "P" stands for "path ordering": for simplicity, we take 5(0, a;) 
along the straight-line path from to a;. 

These field-strength correlators play an important role in hadron physics. In the 
spectrum of heavy QQ bound states, they govern the effect of the gluon condensate on the 
level splittings [p], |2], |3] . They are the basic quantities in models of stochastic confinement 
of colour ||, |5|, H and in the description of high-energy hadron scattering JF|, ||, |], |i0fl . 
In some recent works [JTT , 12| , these correlators have been semi-classically evaluated in 



the single-instanton approximation and in the instanton dilute-gas model, so providing 
useful information about the role of the semiclassical modes forming the QCD vacuum. 

In the Euclidean theory, translational, 0(4)- and parity invariance require the corre- 
lator (1.1) to be of the following form ||, || ||: 

T> mvcT {x) = (dpvdpa - <W) T>(x 2 ) + T>i(x z 

, . dV 1 (x 2 ) 
~\~\x pX v o p a x pXfyd p V -\- x pX a o py x pXyd p a j Qx^ ' V / 

where T> and T>\ are invariant functions of x 2 . 

These functions T>(x 2 ) and T>i(x 2 ) have been directly determined (in the Euclidean 
theory) by numerical simulations on a lattice in the quenched (i.e., pure gauge) theory, 
with gauge-group SU(2) [O, in the quenched SU(3) theory in the range of physical 



distances between 0.1 and 1 fm |14[ [15[ and also in full QCD, i.e., including the effects 



of dynamical fermions [1(J. In another approach fl7f , they have been extracted (in the 



2 



quenched SU(3) theory) from lattice calculations of the heavy-quark potential, under the 
assumptions of the model of the stochastic vacuum j|, || ||. 

It is convenient to define a V\\(x 2 ) and a T>±(x 2 ) as follows: 

V\\ = V + V!+x 2 — 4 , 

ox z 

V± = V + Vx . (1.4) 

In Figs. 1-2 we display the results for T>\\/A\ and T>±/A\ versus the physical distance in 
fermi units, obtained in the quenched (i.e., pure-gauge) theory, with gauge group SU(3): 
data are taken from Refs. fl4|, jlSj . is the fundamental constant of QCD in the 
lattice renormalization scheme, in the pure-gauge case: its value, extracted from the 
string tension fl8| , turns out to be about 4.92 MeV for the gauge group SU(3) and in the 



absence of quarks. 

In Fig. 3 we display the results for T>n/A F and T>±/ ' A F versus the physical distance 
in fermi units, obtained in full QCD, with gauge group SU(3) and Nf — 4 flavours of 
staggered fermions with a quark mass m q = 0.01 in lattice units: data are taken from Ref. 
16]. In Fig. 4 we display the corresponding results for the quark mass m q = 0.02 [16 



Ap is an effective A-parameter for QCD in the lattice renormalization scheme, for the 
gauge group SU(3) and Nf = 4 flavours of quarks. It was defined in Ref. |TB[, where the 
value Ap ~ 1.07 MeV was derived and used. 



In Refs. [H, |I^, |TH| best fits to the lattice data with the functions 



V(x 2 ) = Aoe-^M- + p-e-^ 



\x 



4 



V x {x 2 ) = A ie ~ WXA + ^re- WXa , (1.5) 



x 



were performed: the results obtained for all the various cases are reported in Table I. 
From these results one sees that, in order to obtain a fit with a good x 2 '/N , it is necessary 
to include a perturbative-like term behaving as l/|a;| 4 (indeed, a term of this form is pre- 
dicted by ordinary perturbation theory; see for example Ref. |19| and references therein) 
in addition to the nonperturbative exponential term in the parametrization for both T> 
and V>\. In Table I we also report the results from the fits to the restricted set of data of 
the correlators obtained in Ref. fH| , corresponding to physical distances \x\ > 0.4 fm: in 
this case a fit with only the nonperturbative exponential terms for T> and V>\ turns out 
to be acceptable. 
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In this paper we report the results obtained by fitting the lattice data of the gauge- 
invariant field strength correlators in QCD, displayed in Figs. 1-4, with some particular 
functions which are commonly used in the literature in some phenomenological approaches 
to high-energy hadron-hadron scattering. Therefore, we think that these results may be 
useful to those people working with this alternative parametrization. In Sect. 2 we give a 
brief review on the new parametrization used for the best fits and the results so obtained 
are discussed. In Sect. 3 some quantities of physical interest are extracted from our best 
fits and compared with the corresponding results obtained in the original fits to the lattice 
data. 



2. New fits to the lattice data 



In this section we discuss the results obtained by fitting the lattice data reported in Figs. 
1-4 with the following functions: 

V(x 2 ) = ^kG 2 V{x 2 ) + j^e-W/*- , 

V^x 2 ) = il-^GA^ + ^e-l*, (2.1) 
where V(x 2 ) and VAx 2 ) are so defined: 



2 , / 37r|x| \ „ / 3tt|x| 



2V) = \-£r\Ki 



■sa \ oa 




V^) = {^)K 1 ( 2 m . (2.2) 



3d / \ oa 



Kq and K\ are the modified Bessel functions. These expressions for the "nonperturbative" 
parts of T> and T>i, i.e, the expressions reported in Eqs. (2.1) and (2.2), with the exclusion 
of the "perturbative-like" 1/ 1 a: | 4 pieces, are extensively used in the literature in many 
phenomenological approaches to high-energy hadron-hadron scattering (see for example 
Refs. [0, |20|, |2lJ and references therein). They were proposed for the first time in Ref. 
|T0|| , where also a preliminary fit to lattice data was performed. However, at that time 



only lattice data obtained in the quenched SU(3) theory, in the range of physical distances 
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between 0.4 and 1 fm, were available [14]. Therefore, we think that a re-analysis of this 



parametrization, fitting also the new lattice data now available in the quenched SU(3) 



theory in the range of physical distances between 0.1 and 1 fm |T5| and in full QCD fi~6f , 
may be useful for the practitioners in this field. 

Before discussing the results of the fits, we shall remind the reader of some technical 
details about the parametrization (2.1). The functions T>(x 2 ) and Vi(x 2 ) are normalized 
to 1 in x = 0: 

X>(0) = Pi(0) = 1 . (2.3) 

As we shall discuss in the next section, this implies that the parameter G 2 in Eq. (2.1) 
should be identified with the gluon condensate. The parameter k measures the non- 
Abelian character of the correlator: in fact, k = in an Abelian theory, if there are no 
magnetic monopoles, while there is no reason for the P-term to vanish in a non-Abelian 
theory. The expression for T>(x 2 ) in Eq. (2.2) comes from the following ansatz: 



97 r k 2 



V{x 2 ) = =rra- 2 / d A ke ikx 
64 



91 4 ' 

* 2 + (£) 



(2.4) 



where the length-scale a enters into this parametrization as a "correlation length" , being 
defined as 

r+oo 

a = / d\x\V{x 2 ) . (2.5) 



o 



The correlation function T>(x 2 ) is negative at large distances, with the following asymp- 
totic behaviour: 



The function T>\{x 2 ) is chosen such that 

d 



A + x ll ^-)V 1 (x 2 )=AV(x 2 ) , (2.7) 



dx ilj 



which leads to 



1 



Vi(x 2 ) = T - I L dz 2 V(z 2 ), (2.8) 



and then to the expression in Eq. (2.2). 
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The results of the fits to the lattice data using the parametrization (2.1)-(2.2) for T> 
and T>\ are reported in Table II. The continuum lines in Figs. 1-4 have been obtained 
using the parameters of these best fits, in the cases where all the parameters were free. 
The dashed lines in Figs. 1-2 correspond to the nonperturbative parts only in Eq. (2.1) 
and are derived using the same parameters used for the corresponding continuum lines. 
They are drawn in this particular case {quenched SU(3) theory) in order to illustrate the 
role of the perturbative-like terms. 

From the results in Table II, one sees that, in order to obtain fits with a good x 2 /N, 
it is necessary to include the perturbative-like terms ~ 1/|#| 4 in the parametrization for 
D and T>\. In fact one also sees directly from Figs. 1-2 that these terms are necessary to 
well describe the behaviour of the correlators at small distances (down to 0.1 fm), while 
they are less important in the range of distances between 0.4 and 1 fm. The coefficients 
a and ai turns out to be comparable with (even if slightly smaller than) the coefficients 
derived in the corresponding cases in Table I, obtained using the parametrization (1.5). 
Even restricting the set of quenched data to those obtained in Ref. |14|], corresponding 
to physical distances |a;| > 0.4 fm, one finds that a fit with only the nonperturbative 
terms in Eq. (2.1), i.e., fixing = 0, is more acceptable (when compared to the 

same fit applied to the entire set of data between 0.1 and 1 fm), but the x 2 /N is still 
too high (~ 3.8). When all the parameters in the fit are free, the x 2 /N turns out to be 
acceptable in all the various cases examined, even if it is systematically a bit larger than 
the x 2 /N obtained in the corresponding cases examined in Table I, using the original 
parametrization (1.5). 

Therefore, we can conclude that the expressions (2.1)-(2.2) are a good parametrization 
of the correlators D and T>i, in the range of physical distances where the lattice data are 
available (i.e., 0.1-1 fm for quenched QCD and 0.3-1 fm for full QCD). However, the 
parametrization (1.5) appears to be slightly preferable. In the next section we shall 
discuss some quantities of physical interest which can be extracted from the results of the 
best fits in Tables I and II. 



3. Discussion 

Three quantities of physical interest can be extracted from our fits to the lattice data: 
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(1) The correlation length Iq of the gluon field strengths, defined as ("np" stands for 
"nonperturbative", in the sense explained in the previous sections) 



I)H(O) 



+00 



d\x\V {np \x 2 



(3.1) 



(2) The so-called "gluon condensate" , defined as 



.at 



7T 



47T' 



(3.2) 



(3) The parameter k, defined as 



K 



X>("p)(0) + £>j np) (0) 



(3.3) 



The results obtained are summarized in Table III. 

The quantities Iq and G2 play an important role in phenomenology. The correlation 
length is relevant for the description of vacuum models [|], [5], |J. The relevance of the 
gluon condensate was first pointed out by Shifman, Vainshtein and Zakharov (SVZ) P2"| . 
It is a fundamental quantity in QCD, in the context of the SVZ sum rules. 

The physical meaning of the k parameter has been already discussed in the previous 
section: it measures the non-Abelian character of the correlator, since one expects that 
ft = in an Abelian theory with no magnetic monopoles present. This parameter ap- 
pears explicitely in the parametrization (2.1), by virtue of Eq. (2.3). Instead, using the 
parametrization (1.5) one finds the following expression for k: 



K 



A 



A Q + A ± 



(3.4) 



From the results reported in Table III, one sees that, in both parametrizations (1.5) and 
(2.1)-(2.2), the parameter k appears to decrease when increasing the quark mass, tending 
towards the pure-gauge value (obviously, when evaluating the field strength correlators 
(1.1), the quenched, i.e., pure-gauge, limit coincides with the large quark-mass limit, 
m q — ► 00). In other words, the non-Abelian character of the correlator T> appears to 
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increse when approaching the chiral limit m g — > 0. The values of k obtained using the 
parametrization (2.1)-(2.2) in all the cases examined are a bit larger than the correspond- 
ing values of k obtained using the parametrization (1.5). 

Now let us discuss the results for the correlation length Iq. Using the parametrization 
(1.5) one easily finds Iq = \a, while in the parametrization (2.1)-(2.2) one has Iq = a, 
according to Eqs. (2.5) and (2.3). From the results in Table III, one sees that both A^ 
and a decrease by increasing the quark mass, when going from chiral to quenched QCD. 
Obviously, the difference between A^ and a is due to the different parametrization used 
for the correlators. 



Now we come to the gluon condensate. As pointed out in Ref. the lattice provides 



us with a regularized determination of the correlators. We shall briefly repeat here the 
argumentation originally reported in Ref. flTB|, for the benefit of the reader. At small 



distances x a Wilson operator-product-expansion (OPE) |23| is expected to hold. The 
regularized correlators will then mix to the identity operator 1, to the renormalized local 
operators of dimension four, g f'-G a G a : and mf.q~fqf. (/ = l,...,Nf, Nf being the 
number of quark flavours), and to operators of higher dimension: 

V„{x) ~ Ci{x){l)+C a {x)G 2 + £c f (x)m f {:q f q f :) + ... . (3.5) 



2tt 2 ^ y 's-o ; _, 

The mixing to the identity operator Ci(x) shows up as a c/|x| 4 behaviour at small x. The 
mixings to the operators of dimension four C g (x) and Cf(x) are expected to behave as 
constants for x — > 0, while the other Wilson coefficients in the OPE (3.5) are expected to 
vanish when x — > (for dimensional reasons) . The coefficients of the Wilson expansion are 
usually determined in perturbation theory and are known to be plagued by the so-called 



'infrared renormalons" (see for example Ref. |24j] and references therein). In the same 



spirit of Ref. |16| |, we shall assume that the renormalon ambiguity can be safely neglected 



in the extrapolation for x — > of our correlators. With the normalization of Eq. (3.5), 
this gives C g (0) ~ 1. On the same line, the contribution from the quark operators in (3.5) 



can be neglected, for the reasons explained in Ref. |T(|. Within these approximations, 
one immediately recognizes that the parameter G 2 in the parametrization (2.1)-(2.2), 
with the normalization condition (2.3), coincides with the gluon condensate as defined in 
(3.2). Moreover, when using the parametrization (1.5) for the correlators, one obtains the 



8 



following expression for the gluon condensate: 

G 2 - ^(Ao + A,) . (3.6) 

For both parametrizations (1.5) and (2.1)-(2.2), the gluon condensate G2 appears to in- 
crease with the quark mass, as expected |25| , tending towards the asymptotic (pure-gauge) 
value. However, the values of G2 extracted from the fits using the parametrization (2.1)- 
(2.2) are more than a factor two smaller than the corresponding values of G2 extracted 
from the fits using the parametrization (1.5), in all the cases examined. 
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TABLE CAPTIONS 



Tab. I. Results obtained from a best fit to the data of the gauge-invariant field strength 
correlators with the functions (1.5), in the various cases that we have examined: 
"q" stands for "quenched" data in the range of physical distance between 0.1 and 1 
fm |fT5| ; "q*" stands for "quenched" data in the range of physical distance between 
0.4 and 1 fm 0; "f (I)" and "f (II)" stand for "full-QCD" data with quark mass 
(in lattice units) 0.01 and 0.02 respectively |16|]. An asterisk (*) near the value of 
some parameter means that the parameter was fixed to that value. 

Tab. II. Results obtained from a best fit to the data of the gauge-invariant field strength cor- 
relators with the functions (2.1)-(2.2), in the various cases that we have examined: 
the notation used is the same as in Table I. 

Tab. III. The values of some quantities of physical interest extracted from the best-fit results 
in Tables I ( "exp" ) and II ( "bessel" ) in the cases where all the parameters were free. 
The notation used is the same as in Tables I and II. Reported errors refer only to 
our determination and do not include the uncertainty on the physical scale. 
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Table I 



theory 


A /A 4 
xlO" 8 


^i/A 4 
xlO" 8 


1/(AaA) 


a 




l/(A a A) 


X 2 /N 


q 


3.34(20) 


0.70(10) 


182(3) 


0.69(6) 


0.46(3) 


94(15) 


1.7 


q 


8.37(22) 


2.94(9) 


233(2) 


0(*) 


(*) 


(*) 


33 


q* 


3.11(61) 


0.83(23) 


183(8) 


0.22(12) 


0.12(7) 


0(343) 


1.4 


q* 


3.62(19) 


1.23(7) 


183(3) 


0(*) 


(*) 


(*) 


1.3 


f (I) 


174(24) 


20(10) 


544(27) 


0.71(3) 


0.45(3) 


42(11) 


0.5 


f (I) 


438(17) 


303(9) 


642(8) 


0(*) 


(*) 


0(*) 


51 


f (II) 


348(42) 


46(21) 


631(23) 


0.66(3) 


0.39(3) 


61(20) 


0.7 


f (II) 


734(21) 


354(10) 


713(7) 


0(*) 


0(*) 


0(*) 


27 



Table II 



theory 


G 2 /A 4 
xlO" 8 


K 


l/(oA) 


a 


ai 


l/(A tt A) 


X 2 /N 


q 


0.95(5) 


0.89(2) 


122(2) 


0.56(4) 


0.38(3) 


48(12) 


2 


q 


2.36(5) 


0.82(1) 


150(1) 


(*) 


(*) 


(*) 


54 


q* 


0.62(8) 


0.85(4) 


116(3) 


0.50^ 


n 99 +u.iu 

u - zz -0.04 


0(33) 


1.5 


q* 


1.14(4) 


0.80(1) 


122(1) 


0(*) 


0(*) 


0(*) 


3.8 


f (I) 


53(7) 


0.95(5) 


387(16) 


0.66(2) 


0.41(2) 


0.0(1) 


0.5 


f (I) 


185(4) 


0.60(1) 


457(4) 


(*) 


(*) 


(*) 


84 


f (II) 


97(10) 


0.94(3) 


436(13) 


0.62(3) 


0.36(2) 


0.00(5) 


1 


f (II) 


284(5) 


0.73(1) 


502(3) 


0(*) 


0(*) 


0(*) 


53 



Table III 



theory 


fit 


G 2 
(GeV 4 ) 


K 


Xa, a 
(fm) 


f (I) 


exp 


0.015(3) 


0.90(6) 


0.34(2) 


f (I) 


bessel 


0.007(1) 


0.95(5) 


0.48(2) 


f (II) 


exp 


0.031(5) 


0.88(6) 


0.29(1) 


f (II) 


bessel 


0.0127(13) 


0.94(3) 


0.42(1) 


q 


exp 


0.144(11) 


0.83(3) 


0.220(4) 


q 


bessel 


0.056(3) 


0.89(2) 


0.328(5) 
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FIGURE CAPTIONS 



Fig. 1. The function T>\\/A\ versus the physical distance in fermi units, obtained in the 
quenched (i.e., pure-gauge) theory, with gauge group SU(3) (A L ~ 4.92 MeV). Data 
are taken from Refs. |14], [15), while the curves are obtained from our best fit with 
the functions (2.1)-(2.2) [fit no. 1 in Table II]: the continuum line corresponds to 
the entire correlator, while the dashed line corresponds to its nonperturbative part 
only. 

Fig. 2. The same as in Fig. 1 for the function T>±/A\. 

Fig. 3. The functions T>±/Ap (upper curve) and T>\\/Ap (lower curve) versus the physical 
distance in fermi units, obtained in full QCD, with gauge group SU(3) and Nf = 4 
flavours of staggered fermions with a quark mass m q = 0.01 in lattice units (Ap ~ 
1.07 MeV). Data are taken from Ref. fIE|| , while the curves are obtained from our 
best fit with the functions (2.1)-(2.2) [fit no. 5 in Table II]. 

Fig. 4. The same as in Fig. 3 for a quark mass m q = 0.02 in lattice units (A^ ~ 1.07 MeV). 

Data are taken from Ref. |16[ , while the curves are obtained from our best fit with 
the functions (2.1)-(2.2) [fit no. 7 in Table II]. 
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